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Abstract: We study a large class of reversible Markov chains with discrete state space and 
transition matrix Pj\[. We define the notion of a set of metastable points as a subset of the 
state space such that (i) this set is reached from any point x £ Ftv without return to 
X with probability at least b^, while (ii) for any two point x,y in the metastable set, the 
probability T~y to reach y from x without return to x is smaller than <C ^tv- Under 
some additional non-degeneracy assumption, we show that in such a situation: 

(i) To each metastable point corresponds a metastable state, whose mean exit time can be 
computed precisely. 

(ii) To each metastable point corresponds one simple eigenvalue of 1 — which is essentially 
equal to the inverse mean exit time from this state. Moreover, these results imply very 
sharp uniform control of the deviation of the probability distribution of metastable exit 
times from the exponential distribution. 
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Section 1 



1. Introduction 

In a recent paper [BEGK] we have presented rather sharp estimates on metastable tran- 
sition times, both on the level of their mean values, their Laplace transforms, and their dis- 
tribution, for a class of reversible Markov chains that may best be characterized as random 
walks in multi-well potentials, and that arise naturally in the context of Glauber dynamics 
for certain mean field models. These results allow for a very precise control of the behaviour 
of such processes over very long times. 

In the present paper wc continue our investigation of mctastability in Markov chains fo- 
cusing however on the connection between metastability and spectral theory while working 
in a more general abstract context. Relating metastability to spectral characteristics of the 
Markov generator or transition matrix is in fact a rather old topic. First mathematical results 
go back at least as far as Wentzell [W] and Frcidlin and Wentzell [FW] . Freidlin and Wcntzell 
relate the eigenvalues of the transition matrix of Markov processes with exponentially 
small transition probabilities to exit times from "cycles" ; Wentzell has a similar result for the 
spectral gap in the case of certain diffusion processes. All these relations are on the level of 
logarithmic equivalence, i.e. of the form liniejo ^^^(K'^i) ~ ^ where e is the small parameter, 
and XI, are the eigenvalues, resp. exit times. For more recent results of this type, see 
[M,Sc]. Rather recently, Gaveau and Schulman [GS] (see also [BK] for an interesting discus- 
sion) have developed a more general program to give a spectral definition of metastability in 
a rather general setting of Markov chains with discrete state space. In their approach low 
lying eigenvalues are related to metastable time scales and the corresponding eigenfunctions 
are related to metastable states. This interesting approach still suffers, however, from rather 
imprecise relations between eigenvalues and time-scales, and eigenfunctions and states. 

In this paper wc will put these notions on a mathematically clean and precise basis for 
a wide class of Markov chains Xf with countable state space r^r^, indexed by some large 
parameter N. Our starting point will be the definition of a metastable set of points each 
of which is supposed to be a representative of one metastable state, on a chosen time scale. 
It is important that our approach allows to consider the case where the cardinality of A^at 
depends on N. The key idea behind our definition will be that it ensures that the time it 
takes to visit the representative point once the process enters a "metastable state" is very 
short compared to the lifetime of the metastable state. Thus, observing the visits of the 



We expect that this approach can be extended with suitable modifications to processes with continuous 
state space. Work on this problem is in progress. 
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process at the metastable set suffices largely to trace the history of the process. We will then 
show that (under certain conditions ensuring the simplicity of the low-lying spectrum) the 
expected times of transitions from each such metastable point to "more stable" ones (this 
notion will be defined precisely later) are precisely equal to the inverse of one eigenvalue (i.e. 
Ti = A~^(l + 0(1))) and that the corresponding cigcnfunction is essentially the indicator 
function of the attractor of the corresponding metastable point. This relation between times 
and eigenvalues can be considered as the analogue of a quantum mechanical "uncertainty 
principle". Moreover, we will give precise formulas expressing these metastable transition 
times in terms of escape probabilities and the invariant measure. Finally, we will derive 
uniform convergence results for the probability distribution of these times to the exponential 
distribution. Let us note that one main clue to the precise uncertainty principle is that we 
consider transition times between metastable points, rather than exit times from domains. 
In the existing literature, the problem of transitions between states involving the passage 
through some "saddle point" (or "bottle neck") is almost persistently avoided (for reasons 
that we have pointed out in the introduction of [BEGK]), except in one-dimensional situations 
where special methods can be used (as mentioned e.g. in the very recent paper [GM]). But 
the passage through the saddle point has a significant impact on the transition time which 
in general can be neglected only on the level of logarithmic equivalence^. Our results here, 
together with those in [BEGK], appear to be the first that systematically control these effects. 

Let us now introduce our setting. We consider a discrete time^ and specify our Markov 
chains by their transition matrix P/v whose elements p]\r{x,y), a;,y € Fat denote the one-step 
transition probabilities of the chain. In this paper we focus on the case where the chain is 
reversible^ with respect to some probability measure Qat on Fjy. We will always be interested 
in the case where the cardinality of Fjv is finite but tends to infinity as | cx). Intuitively, 
metastability corresponds to a situation where the state space F^v can be decomposed into 
a number of disjoint components each containing a state such that the time to reach one of 
these states from anywhere is much smaller than the time it takes to travel between any two 
of these states. We will now make this notion precise. Recall from [BEGK] the notation rf 
for the first instance the chain starting in x at time reaches the set / C Fjv, 

Tf =mi{t>0:Xte l\Xo = x} (1.1) 



E.g. the lack of precision in the relation Tm = C(l/(1 — (1 — A)*)) in [GS] is partly due to this fact. 
There is no difficulty in applying our results to continuous time chains by using suitable embeddings. 
'The case of irreversible Markov chains will be studied in a forthcoming publication [EK]. 
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DeRnition 1.1: A set A^at C Ftv will be called a set 0/ met ast able points if it satisfies 
the following assumptions. For finite positive constants qn, such that, for some sequence 
s:n i 0, a'^ < e^bN it holds that 



(i) For allzer 



N, 



(a) For any x y E Mn, 



[r^, < rl] > bN (1.2) 



[t^ < r:] < a-' (1.3) 



We associate with each x G A^at its local valley 



Aix) ^ G r^v : P = rj,^] = sup P [r^ = rj,^] | (1.4) 



We win set 



and 



- (1-5) 

rjv = max Rx <^ 

_i (1-6) 
c»r = min i?^. > 

Note that the sets A{x) are not necessarily disjoint. We will however show later that the 
set of points that belong to more than one local valley has very small mass under Qjy. The 
above conditions do not fix at uniquely. It will be reasonable to choose M. jv always such 
that for all X G AIat, 

QAr(x)= sup QAr(z) (1.7) 

zeA(a:) 

The quantities P [rf < r^], 1 C A^at furnish crucial characteristics of the chain. We will 
therefore introduce some special notation for them: for I C TWa^ and x G M.n\I, set 

Tx,i = (P[rf < T^])-' (1.8) 

and 

Ti= sup Txj (1.9) 

xeMN\i 

Note that these quantities depend on N, even though this is suppressed in the notation. 
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For simplicity we will consider in this paper only chains that satisfy an additional assump- 
tion of non- degeneracy: 

Definition 1.2: We say that the family of Markov chains is generic on the level of the set 
Mn, if there exists a sequence e^r i 0, such that 

(i) For all pairs x,y & Mn, o,nd any set I C M.N\{x,y} either T^j < e^Tyj or Tyj < 

(a) There exists mi G Mn, s.t. for all x G MN\mi, Q.n{x) < ejvQAr(Tni). 

We can now state our main results. We do this in a slightly simplified form; more precise 
statements, containing explicit estimates of the error terms, will be formulated in the later 
sections. 

Theorem 1.3: Consider a discrete time Markov chain with state space Fn, transition 
matrix Pn, and metastable set Mn (as defined in Definition 1.1). Assume that the chain is 
generic on the IcvcIMn in the sense of Definition 2.1. Assume further thatrN£N\F n\\Mn\ i 
0, and rj^CN^N i 0, as N oo. For every x E Mn set Mn{x) = {y G Mn '■ Qjv(y) > 
Qn{x)}, define the metastable exit time tx = ^XfivCa:)' ^'^^^ 

(i) For any x E Mn, 

Etx = R-^Tx,M^(^.)i'^ + o{l)) (1.10) 

(a) For any x E Mn, there exists an eigenvalue Xx of 1 — Pn that satisfies 

A. = ^(1 + 0(1)) (1.11) 
Moreover, there exists a constant c > such that for all N 

a{l - Pn)\ IJxeMN K C (c6jv|rjvr\ 1] (1.12) 
(here a{l — Pn) denotes the spectrum of 1 — Pn )■ 

(Hi) If 4>x denotes the right- eigenvector of Pn corresponding to the eigenvalue Xx, normalized 
so that 4>x{x) = 1, then 

1 O(ejv), otherwise 

(iv) For any x G Mn, for all t > 0, 

¥[tx > tEtx] = e-*^^+°^^'>\l + o{l)) (1.14) 
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Remark: We will see that P[t| < t^^^^^] is extremely close to one for all y G A{x), with 
the possible exception of some points for which Qjv(y) Qn{x)- Therefore, the correspond- 
ing (normalized) left eigenvectors ipx {y) = "^^^^^^''^^^f"! are to very good approximation 
equal to the invariant measure conditioned on the valley A{x). As the invariant measure Qn 
conditioned on A{x) can be reasonably identified with a metastable state, this establishes 
in a precise way the relation between eigenvectors and metastable distributions. Brought 
to a point, our theorem then says that the left eigenfunctions of 1 — P/v are the metastable 
states, the corresponding eigenvalues the mean lifetime of these states which can be com- 
puted in terms of exit probabilities via (1.10), and that the lifetime of a metastable state is 
exponentially distributed. 

Remark: Theorem 1.3 actually holds under slightly weaker hypothesis than those stated 
in Definition 1.2. Namely, as will become clear in the proof given in Section 5, the non- 
degeneracy of the quantities T^j is needed only for certain sets /. On the other hand, if these 
weaker conditions fail, the theorem will no longer be true in this simple form. Namely, in a 
situation where certain subsets Si C Mn are such that for all x G Si, T^j (for certain relevant 
sets /, sec Section 5) differ only by constant factors, the eigenvalues and eigenfunctions 
corresponding to this set will have to be computed specially through a finite dimensional, 
non-trivial diagonalisation problem. While this can in principle be done on the basis of the 
methods presented here, we prefer to stay within the context of the more transparent generic 
situation for the purposes of this paper. Even more interesting situations crating genuinely 
new effect occur when degenerate subsets of states whose cardinality tends to infinity with 
N are present. While these fall beyond the scope of the present paper, the tools provided 
here and in [BEGK] can still of use, as is shown in [BEG]. 

Let us comment on the general motivation behind the formulation of Theorem 1.3. The 
theorem allows, in a very general setting, to reduce all relevant quantities governing the 
metastable behaviour of a Markov chain to the computation of the key parameters, T^^y and 
Rxi x,y G A4n- The first point to observe is that these quantities are in many situations 
rather easy to control with good precision. In fact, control of requires only knowledge 
of the invariant measure. Moreover, the "escape probabilities", T~y, arc related by a factor 
QAr(x) to the Newtonian capacity of the point y relative to x and thus satisfy a variational 
principle that allows to express them in terms of certain constraint minima of the Dirichlet 
form of the Markov chain in question. In [BEGK] we have shown how this well-known fact 
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(see e.g. [Li], Section 6) can be used to give very sharp estimates on these quantities for 
the discrete diffusion processes studied there. Similar ideas may be used in a wide variety of 
situations (for another example, see [BBG]); we remind the reader that the same variational 
representation is at the basis of the "electric network" method [BS]. Let us mention that 
our general obsession with sharp results is motivated mainly by applications to disordered 
models there the transition matrix Pjv is itself a random variable. Fluctuation effects on the 
long-time behaviour provoked by the disorder can then only be analysed if sharp estimates 
on the relevant quantities are available. For examples see [BEGK, BBG]. 

In fact, in the setting of [BEGK], i.e. a random walk on {Z/N)'^nA with reversible measure 
Qn{x) = exp(— iVFiv(a;)), where Fn is "close" to some smooth function F with finite number 
of local minima satisfying some additional genericity requirements, and the natural choice for 
Mn being the set of local minima of Fn, the key quantities of Theorem L3 were estimated 
as 

bN > cN-^/'^ (L15) 
rjv < cN-'^/^ 

(L16) 

CN < CN"/^ 

T^y = e^(i) jv-(''-2)/2e^[^'v(^*(^'2'))--P'^W] (1.17) 

where z*{x,y) is the position of the saddle point between x and y. Moreover, under the 
genericity assumption of [BEGK], 

Civ < e"^" (1.18) 

for some a > 0. The reader will check that Theorem 1.3, together with the precisions detailed 
in the later sections, provides very sharp estimates on the low-lying eigenvalues of 1 — and 
considerably sharpens the estimates on the distribution function of the metastable transition 
times given in [BEGK]. 

Let us note that Theorem 1.3 allows to get results under much milder regularity assump- 
tions on the functions Fn then were assumed in [BEGK]; in particular, it is clear that one 
can deal with situations where an unbounded number of "shallow" local minima is present. 
Most of such minima can simply be ignored in the definition of the metastable set M n which 
then will take into account only sufficiently deep minima. This is an important point in many 
applications, e.g. to spin glass like models (but also molecular dynamics, as discussed below), 
where the number of local minima is expected to be very large (e.g. exp(oA'')), while the 
metastable behaviour is dominated by much fewer "valleys" . For a discussion from a physics 
point of view, see e.g. [BK]. 
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A second motivation for Theorem 1.3 is given by recent work of Schiitte et al. [S,SFHD]. 
There, a numerical method for the analysis of metastablc conformational states of macro- 
molecules is proposed that relies on the numerical investigation of the Gibbs distribution 
for the molecular equilibrium state via a Markovian molccTilar dynamics (on a discretized 
state space) . The key idea of the approach is to replace the time-consuming full simulation 
of the chain by a numerical computation of the low-lying spectrum and the corresponding 
eigenf unctions, and to deduce from here results on the metastable states and their life times. 
Our theorem allows to rigorously justify these deductions in a quantitative way in a setting 
that is sufficiently general to incorporate their situations. 

The remainder of this article is organized as follows. In Section 2 we recall some basic 
notions, and more importantly, show that the knowledge of T^^y for all x,y e Mn is enough 
to estimate more general transition probabilities. As a byproduct, we will show the existence 
of a natural "valley-structure" on the state space, and the existence of a natural (asymptotic) 
ultra-metric on the set Mn- In Section 3 we show how to estimate mean transition times. 
The key result will be Theorem 3.5 which will imply the first assertion of Theorem 1.3. In 
Section 4 we begin our investigation of the relation between spectra and transition times. 
The key result there is a characterization of parts of the spectrum of (1 — P/v) in terms 
of the roots of some non-linear equation involving certain Laplace transforms of transitions 
times, as well as a representation of the corresponding eigenvectors in terms of such Laplace 
transforms. This together with some analysis of the properties of these Laplace transforms 
and an upper bound, using a Donsker-Varadhan [DV] argument, will give sharp two-sided 
estimates on the first eigenvalue of general Dirichlct operators in terms of mean exit times. 
These estimates will furnish a crucial input for Section 5 where wc will prove that the low-lying 
eigenvalues of 1 — Pjy are very close to the principal eigenvalues of certain Dirichlct operators 
(1 — Pn)^^ , with suitably constructed exclusion sets Sj. This will prove the second assertion 
of Theorem 1.3. In the course of the proof wc will also provide rather precise estimates on 
the corresponding eigenfunction. In the last Section wc use the spectral information obtained 
before to derive, using Laplace inversion formulas, very sharp estimates on the probability 
distributions of transition times. These will in particular imply the last assertion of Theorem 
1.3. 

Acknowledgements: We would like to thank Christof Schiitte and his collaborators for 
explaining their approach to conformational dynamics and very motivating discussions. 
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2. Some notation and elementary facts. 

In this section we collect some useful notations and a number of more or less simple facts 
that we will come back to repeatedly. 

The most common notion we will use are the stopping times rf defined in (1.1). To avoid 
having to distinguish cases where x G I, it will sometimes be convenient to use the alternative 
quantities 

af = mm{t > : Xt e I\Xo = x] (2.1) 
that take the value if x G /. 

Our analysis is largely based on the study of Laplace transforms of transition times. For 
IcTn we denote by (Pn)^ the Dirichlet operator 

{PnY = I/cPjv : I/cM^^ I/cM^'^, r = rjv\i" (2.2) 

Since our Markov chains are reversible with respect to the measure Qn, the matrix (Pn)^ is 
a symmetric operator on I/c^^(rjv, Qat) and thus 

||(Pjv)^|| = max{|A| I A G (tUPnY)} (2.3) 

where 1 1 • || denotes the operator norm induced by ]l7c^^(rjv, Qat)- For a point x G Fjv, subsets 
/, JcFjv and u G C, ^{u) < - log ||(Pjv)^^''||, we define 

OO 

GIj{u) ^E[e--' = ^e-*F[rf = t < r]] (2.4) 

t=i 

and 

( GJj{u) for x^lUJ, 
Kf^j{u)=E[e'"''l^.<^.] = I 'l for x £ I, (2.5) 

( for X G J\/ 

The Perron-Frobenius theorem applied to the positive matrix (PnY implies that G^ j{u) and 
K'f J {u) converge locally uniformly on their domain of definition, more precisely 

- log I |(P7v)^| I = svi^{u G M I i^7,7(w) exists for all x i I] (2.6) 

We now collect a number of useful standard results that follow trivially from the strong 
Markov property and/or reversibility, for easy reference. 

From the strong Markov property one gets: 
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Lemma 2.1: Fix I, J, LcTn- Then for all ^{u) < - log ||(Pjv)^'-'-^|| 

GlAu) = G?\i,juLN + E G^yJuJuMKlAu), xeV^ (2.7) 

yeL 

In the following we will adopt the (slightly awkward) notation PnF^ = YlzeVN ^n{x, z)F^ 
The following are useful specializations of the foregoing result which we state without proof: 

Corollary 2.2: FixI,JcTN. Then forxeVN 

e^PNKfjiu) = GIj{u), xETn (2.8) 

and 

(1 - e^PN)duKlj{u) = GIj{u), x^IUJ (2.9) 
where 9„ denotes differentiation w.r.t. u. 

The following renewal equation will be used heavily: 

Corollary 2.3: Let IcTn- Then for allx^IUy and Sft(u) < - log IKPat)^*^^!! 

finally, from reversibility of the chain one has 
Lemma 2.4: Fix x,y eTn and IcFn- Then 

Qiv(a^)Q7ux = Qiv(y)G'^,/u, (2-11) 



The next few Lemmata imply the existence of a nested valley structure and that the 
knowledge of the quantities T^^y and the invariant measure are enough to control all transition 
probabilities with sufficient precision. The main result is an approximate ultra-metric triangle 
inequality. Let us define (the capacity of x relative to y) E{x,y) = QN{x)T~y. We will show 
that 

Lemma 2.5: Assume that y,m € and JcT]\i\y\m such that for < 5 < E{m, J) < 
6E(m,y). Then 

1-26 E(m,J) 1 

< — ^ — — < (2.12) 

1-S - E{y,J) -1-5 ^ ^ 
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Proof: We first prove the upper bound. We write 

xeJ ' 

Now 



Fir'" < r™ 

[r^ < r]] = P[< < rl, r,- < r]] + F[t^ < r]^^]J^J—-^ (2.14) 



[rTuy < ^m] 



Now by assumption 

^^''J ^ 'yUm\ . ir-lTJ 



< J < S (2.15) 



Inserting (2.15) into (2.14) we arrive at 

< r^j] < Ar^ <r^j,r^„,< r^] + <5P[r;^ < r^^ J < P[r,- < rj] + 5¥[t^^ < r^j] (2.16) 
Inserting this inequahty into (2.13) imphes 

< C] < (1 - ^r'^ri^nr'j < r^] (2.17) 
We now turn to the lower bound. We first show that the assumption imphes 

nr'j < rl] < 6(1 - S)-' (2.18) 

Namely, 

< C] > <rT< C] = P[r™ < r^^^jPfrJ < r^] (2.19) 

< rTum] =IPkr < C] - < rT < C] 

>P[r™<C]-P[T7<C] (2.20) 

>P[r-<C](l-5) 
where the last inequality follows from the assumption. Thus 

nrT < C] > nrT < < T^](l - S) (2.21) 

Solving this inequality for P[tj < r^], the assumption yields (2.18). 

We continue as in the proof of the upper bound and write for x e J, using (2.18), 
P[r; < t]] = F[r^ < rj, r^, < rj] + P[t,- < tJu^]F[t^ < r^] 

(2 22) 

<P[r;;^<T^]+F[r,-<r^]5(l-5)-i 



But 
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proving 

F[t^ < T^j] < P[r^ < rj] (2.23) 

Inserting (2.23) into (2.13) for m = y and, using once more (2.13) in the resulting estimate, 
we obtain 

which yields the lower bound in (2.12). 
Corollary 2.6: Assume that x,y,z e Mn- Then 

E{x, y)>^ min iE{x, z),E{z, y)) (2.25) 

Proof: By contradiction. Assume that E{x,y) < ^ min (E{x, z), E{z,y)). Then E{x,y) < 
^E{x,z), and so by Lemma 2.5, 

1 E(x,y) 3 

- < ^) [ < - 2.26 
2-E{z,y)-2 ^ ' 

and in particular E{y,z) < 2E{x,y), in contradiction with the assumption. •() 
If we set 

f ~\nE(x,y), if X ^ y , ^ 

e{x,y)^\ ^ (2.27) 

10, li x = y 

then Lemma 2.5 implies that e furnishes an "almost" ultra-metric, i.e. it holds that e{x,y) < 
max(e(a;, z), e{z, y)) + ln3 which will turn out to be a useful tool later. We mention that in 
the case of discrete diffusions in potentials, the quantities e{x, y) are essentially N times the 
heights of the essential saddles between points x and y. 

The appearance of a natural ultra-metric structure on the set of metastable states under 
our minimal assumptions is interesting in itself. 

A simple corollary of Lemma 2.5 shows that the notion of elementary valleys, A{m), is 
reasonable in the sense that "few" points may belong to more than one valley. 

Liemma 2.7: Assume that x,m E Mn and y eTn- Then 

nrl<ry]>e and P[t,^ < r,^] > e (2.28) 

implies that 

Qiv(y) < 2e-iQ;v(m)P[C < C] (2-29) 
We leave the easy proof to the reader. 
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In this chapter we will prove various estimates of conditioned transition times E[rf |t/ < 
Tj], where / U JcMn- The control obtained is crucial for the investigation of the low lying 
spectrum in Chapters 4 and 5. In the particular setting of the paper [BEGK], essentially the 
same types of estimates have been proven. Apart from re-proving these in the more abstract 
setting we consider here, we also present entirely different proofs that avoid the inductive 
structure of the proofs given in [BEGK]. Instead, it uses heavily a representation formula for 
the Green's function (which first appeared in Section 3, Eq. (3.12) of [BEGK]^). While the 
new proofs are maybe less intuitive from a probabilistic point of view, they are considerably 
simpler. 

Theorem 3.1: Fix a nonempty, irreducible, proper subset QcT^- Let (1 — Pn)^ denote 
the Dirichlet operator with zero boundary conditions at Q'^. Then the Green's function defined 
as G%{x,y) = ((1 - PN)^^)~'^'S.y{x), x,y e^, is given by 



Proof: This theorem follows essentially from the proof of Eq. (3.12) of [BEGK]. Using e.g. 
the maximum principle, it follows that (1 — Pn)^" is invertible. Prom (2.8) we obtain, using 
(2.5), 

(1 - P^f^Kl^^O) = My)Gf,.jO) {x,ye U) (3.2) 

This function serves as a fundamental solution and we compute for x,y G Q, using the 
symmetry of (1 — Pn)^" , 

QN{x)Gf,.JO)G%\x,y)={{l-PNf^K^^'l^.{0),G%'{;y))Q, 

= (<},.(0),(1 -P7v)^^G^°(-,y))Q. (3.3) 

This proves (3.1).0 

Remark: Observe that (3.1) still makes sense for x G and y € dil,, where we define the 
boundary dl of a set IcTn to be 

dl = {xer\3yel: PN{y, x) > 0} (3.4) 

^More recently, the same formula was rederived by Gaveau and Moreau [GM] also for the non- reversible 
case. 
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For such X and y reversibility (2.11) and the renewal relation (2.10) for -u = and 7 = 0'^ 
imply 

{x, y) = P[r^- = T^e] (x G 1^, y G d^) (3.5) 

Based on Theorem 3.1 we can derive an alternative representation of a particular h- 
transform of the Green's function with h{y) = P[t|' < Tj] that will prove useful in the 
sequel. 

Proposition 3. 2: For every nontrivial partition JU J = O'^ such that I and J are not empty 
and I\J communicates with O we have 

F[rf < Tj]-^G%\x,y)F[Tf < r^j] = ^i^Ll^L^L^^Anc(x, y), x,y G n (3.6) 



where 



Furthermore, 



Anc(x,y) = —, — 777, x,yEO, (3.7) 



^<Anc(x,y)<3 (3.8) 



Proof: (3.6) is a straightforward calculation that uses the renewal equation (2.10), reversibil- 
ity, and the strong Markov property. Indeed, by (3.1) the left-hand side of (3.6) equals 



Q^(y)P[ag<r^.MrJ<r5] 
Qjv(x)P[a^e < r-]P[r/ < r]] 

By the renewal equation, this equals 

Q^r{y)F[ay<ry,.^^]F[Tf<ry] 



(3.9) 



Qn{x)F[t^.^^ < r,^]P[c75. < T-]F[Tf < T-j] 
which by reversibility turns into 



(3.10) 



F[a- < Tf\Tf < T^lPKeu. < ry]F[Ty < r^.] 



(3.11) 
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where the last identity uses that by the strong Markov property 

< ^f^^i ^ = < ^ = < ^fuAnrf < r]] (3.12) 
(3.11) immediately implies (3.6). 

We now turn to the proof of the bound (3.8). Since Aqo(x, x) = 1 it is enough to consider 
the case where x ^ y. Moreover, since lS.^c(^x,y) = A^e^y x) ' upper bound Anc(x, y) < 3 
will immediately imply the claimed lower bound. 

The basic input here is the observation that a path from y to VL^ either visits a point x or 
it does not, yielding, together with the strong Markov property 

Using this identity for the first factor in the numerator of (3.7), we obtain that Anc(a;, y) can 
be written as Anc(a;, y) = (/) + (//) where 

The renewal equation was used in the second equality. Decompose the event in the second 
factor of the numerator and use (3.13) in the denominator. This yields 

^ ^ " F[ra.<r.^u.]+P[r,^<r^cuJP[r^.<r|] ^ ^^^^^ < ^ J + ^ ^ ^ (3.15) 
For (II) we get 

_ < T^^uyMr^^ < ^xujP[r^^u, < rg] _ P[t| < r^eujP[rf,. < rguj ^ 

^ ^ P[t^cu, < r-]P[T^e < T-]P[r^eu. < Pfe < r-nr^o^x < r^] ' ^' ^ 

The bounds (3.8) are now obvious. <0> 

The representation (3.6) for the Green's function implies immediately a corresponding 
representation for the (conditioned) expectation of entrance times rf . To see this, recall 
from (2.9) for u = that 

[l-P^Y^'^EUl^^y^^yy] = F[t^ < r% y^lUJ (3.17) 



This yields immediately the 
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Corollary 3.3: Let I, JcTn- Then for allx^lUJ 

nrf\Tf<r^j]= nrf<T]]-'G%\x,y)nrf<ry] 



A first consequence of the representation given above is 
Corollary 3.4: Fix IcMn- Then for all x eTn 

m\Tf <T^^^^j]<3b^'\Tj,\ (3.19) 

In particular, 

nrX,J<3b-'\r^\ (3.20) 
Proof: Using (3.7) in (3.19), we get that 

Using the lower bound (1.2) from Definition 1.1 together with the upper bound (3.8), we get 

E [rf\rf < TX,,\r] < 3&]v' E < 1^/ < ^Xl.X/] (3-22) 

yerN\MN 

from which the claimed estimate follows by bounding the conditional probability by one^°. 
The special case I = Mn follows in the same way, with the more explicit bound 

ErX,, <36^' E n^l<rl,,\ (3.23) 
y&TN\MN 

This concludes the proof of the corollary. <(> 

Theorem 2.2 allows to compute very easily the mean times of metastable transitions. 
Theorem 3.5: Assume that J C M.n , x & Ai^ , and x,J satisfy the condition 

T^,j = Tj (3.24) 



^'^It is obvious that in cases when iFjv | = oo this bound can in many cases be improved to yield a reasonable 
estimate. Details will however depend upon assumptions on the global geometry. 
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We now consider the remaining contributions. This is bounded by 

1 



J2 ^rn (3.34) 



where 

L^= L^{y)^ E '^NiynolKT^j] (3.35) 

yeA{m)\J veA{m)\J 

Assume first that y is such that 
(CJ) Qiv(y)P[T^ < ry] ~ Qiv(m)P[r7 < r^] and 
(Cx) QAr(y)P[T| < r^^] ~ Qiv(m)P[T- < r™] hold, 

where we introduced the notation a~644>i<|<3. Then 

L^(y)<9Q^.(y)pL|,^ (3.36) 

There are two cases: 



(i) If E{m, J) < \E{m, x), then by Lemma 2.5, ^Q^jg)pi;i<;i]^ < | or 



Qjv(m) < ^Qjv(x)^ < ejv^Qiv(a;) (3.37) 

Hence 

^m(2/) < Qiv(y) < ^^^eN^R^'^NiAix)) (3.38) 

(ii) liE{m,J) > ^E(m,x), then E{x,J) > lE{m,x) or Q7v(a;)P[Tj < r^] > |QAf(m)P[r^ < 
r™] so that 

< ^Q^MQid^^ ^ „,,n.^MA(,)) (3.39) 

Finally we must consider the cases where (CJ) or (Cx) are violated. 

(iii) Assume that (Cx) fails. Then by Lemma 2.5, P[r^ < t^] > ^Pir™ < r™] which implies 
that 

Lmiy) < Qiv(y) < 3Qiv(m)5rV4^ ^ 3QAr(m)^i^^^^-^^l^ 

On On (In 
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(iv) Finally it remains the case where (CJ) fails but (Cx) holds. Then P[rJ < r-j] > ^Pir^ < 
ry] > ^ and Qjv(y)P[r^ < r^] < fQiv(m)P[r- < C] = IQn{x)F[t?^ < t|]. Thus 
Lmiy) satisfies equally the bound (3.40). 

Using these four bounds, summing over y one gets 



.1 \ MN\\A{m )\ 
6 AT Oat 



Lm < 27qN{A{x)) max ( eNR.R-\ ' ^ [. ) (3-41) 



Putting everything together, we arrive at the assertion of the theorem. <0> 

Remark: As a trivial corollary from the proof of Theorem 3.5 one has 

Corollary 3.6: Let x ^ Mn and J <Z Mn{x). Then the conclusions of Theorem 3.5 also 
hold. 

Finally, we can easily prove a general upper bound on any conditional expectation. 
Theorem 3.7: For any x G Fjv and I,JcMn, 

E [rflrf <T]]<C sup {RmnrZj < CD"' (3-42) 

meMN\I\J 

To prove this theorem the representation of the Green's function given in Proposition 2.2 
is particularly convenient. It yields 

nrf\rf<r]]= ^1 pr^/ ' ^1 ^ ^^^(^'^) (3-43) 

Note first that the terms with y such that P[t]'u j < t^] > Sb^ yield a contribution of no more 
than \T]^\{5bf^)~^ which is negligible. To treat the remaining terms, wc use that whenever 
y € A{m), Lemma 2.5 implies that Pir/^j < r^] > ^^P[r]^j < r™]. Thus 

meMN\i\J yeA{m) ; I Jul mJ ^^^^^ 

SlFivl op-1 1 



from which the claim of the theorem follows by our general assumptions. Note that by very 
much the same arguments as used before, it is possible to prove that 

< ^ ^ (1 + < r^rf < r]] (3.45) 
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which ahows to get the sharper estimate 

E[.f|xf<rJl<f5l+ E Hl + i)R-J?^&±lfK^ (3.46) 



meMN\I\J ^ •^'-^^ 







We conclude this section by stating some consequences of the two preceding theorems that 
win be useful later. 

Lemma 3.8: Let I,m satisfy the hypothesis of Theorem 3.5. Then 

maxE[rJ] = E[rr](l + 0{Tium/Ti)) (3.47) 

E[r-,r™<rf 



Moreover, we have 



E[Tr](l-C>(rjum/r7))) (3.48) 



In particular, 



C < rT\ = R^n' (1 + ©(T/um/T/)) (3.49) 



Proof: Decomposing into the events where m is and is not visited before /, and, using the 
strong Markov property, one gets 

E[rf] = ¥[Tf < r^nrflrf < +P[C < rf](E[C|r^ < r/] +E[rr]) (3.50) 

Using Theorems 3.5 and 3.7, this implies (3.47) readily. In the same way, or by differentiating 
the renewal equation (2.10), one gets 

]pr_m _m _ ml 

E[tj J = E[tj \tj < r^J + p[^m ^ ^m] (^-51) 

Bounding the first summand on the right by Theorem 3.7 gives (3.48). Using Theorem 3.5 
for the right hand side of (3.48) gives (3.49). 
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4. Laplace transforms and spectra 

In this section we present a characterization of the spectrum of the Dirichlet operator 
(1 — PnY , IcMn, in terms of Laplace transforms of transition times (defined in (2.4) and 
(2.5)). This connection forms the basis of the investigation of the low-lying spectrum that is 
presented in Section 5. To exploit this characterization we study the region of analyticity and 
boundedness of Laplace transforms. As a first consequence we then show that the principal 
eigenvalue for Dirichlet operators are with high precision equal to the inverse of expected 
transition times. A combination of these results then leads to the characterization of the 
low-lying spectrum given in the next section. 

For any JcAi n we denote the principal eigenvalue of the Dirichlet-operator by 

Aj = min(7((l-PAr)-^) (4.1) 

For /, JczMn we define the matrix 

GlAn)^(Sm',m-G^j^j{u)) ^ (4.2) 
V / m' ,mEJ\I 

where Sx,y is Kronecker's symbol. We then have 

Lemma 4.1: Fix subsets I, JcMn such that J\I ^ and a number < A = 1 — e~" < 
A/uJ- Then 

A G (t((1 - Piv)0 ^ deigi,j{u)=0 (4.3) 
Moreover, the map kev gi^j{u) 3 (/) ^ (/) E TLjcM^^ defined by 

(j){x)= ^ ^mK^juji^), xeTN (4.4) 

meJ\I 

is an isomorphism onto the eigenspace corresponding to the eigenvalue A. 

Proof: Assume that (p is an eigenfunction with corresponding eigenvalue A < A/yj- We 
have to prove that Qi^j{u) is singular. In view of (2.6) the condition A/uJ > implies that 
defined below is finite. 

4>=Y1 '^("^)^m,/uj(«)> xETn (4.5) 
Furthermore, (2.8) and (2.5) imply for x G Fjv 

e"(l - Pat - (1 - e-"))0(x) = (1 - e"PA,)0(x) = ^ 5w,x ('^-'-^ " ^-!/ujM) 

m'eHJJ meJ 

(4.6) 
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Let A = (p — (j). We want to show A = 0. Obviously, we have A vanishes on / U J and (j) on 
/. Combining (4.6) with the eigenvalue equation for (f) and the choice of u, we obtain 

(1 - Pn^^'A = I(,u j)c (1 - PnYA = I(,u j)c ((1 - PnYcP - (1 - Pn)4>) 
= %uj)c(A(/)-(l-e-")^) = AA 

Since A ^ (7((1— Pjv)^*^"^), we conclude A = 0. Replacing ^ by <^ in (4.6) and, using A = 1— e~" 
again, gives 

0= ^m',.J2^('^^{^rn',m-G^j^j{u)), X^F (4.8) 

m'elUJ meJ 

Choosing x € J\I yields that {(l){m))^^j\^j G kcr Qj^j{u) and the right-hand side of the 
equivalence in (4.3) follows. In particular, we have proven that the restriction map (p i— > 
((/)(m))^gjy defined on the eigenspace corresponding to A is the inverse of the map defined 
in (4.4). 

For the converse implication we note that for A < Xjuj the entries of the matrix Qi,j{u) are 
finite. We replace {<i>{'m))^^j\i in (4.5) by the solution of the linear system Qij{u)(f) = 
and deduce from (4.6) and (4.8) that A is an eigenvalue with eigenfunction 

As a first step we now derive a lower bound on these eigenvalues, using a Donsker-Varadhan 
[DV] like argument that we will later prove to be sharp. 

Lemma. 4.2: For every nonempty subset JgMn we have 

AjmaxE[Tj]>l (4.9) 



Proof: For G M^'v ^^ve for all x, y G Fat and C > 

mm < limfc + Hvf/c) (4.10) 

Thus choosing C = '>p{y)/'4'ix), where tp G MJ"^ is such that ip{x) > for all x G supp^, we 
compute, using reversibility, 

{PNct>,cp)Q^ <l J2 QN{x)PN{x,ymxfm)/m)+m\m/i^{y))) 
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Let be an eigenftmction for the principal eigenvalue and set il^{x) = IE[crj], x G Fjv- Invoking 
(2.9) for u = and / = J we get 

AjIHIL>('^M</')q. (4.12) 

which in turn gives the assertion. •(> 

Wc now study the behavior of Laplace transforms slightly away from their first pole on 
the real axis. 

Lemma 4.3: Fix nonempty subsets I, JcMn- Let Gf j be the Laplace transform defined 
in (2.4)- It follows that for some c > and for A; = 0, 1 uniformly in < 5R(n), < 
c/{cnThjj) and x eTn 

duGlAu) = (1 + 0{\u\cr,Tjuj)) d'^GlAO) (4.13) 



Proof: By (2.6), we know that GJ j{u), x G Tjv, are finite for all u such that 1 — e ^^"^ < 
A/uJ- Put 

Ku,v^Ki-^j{u)-K'i-^j{v) (4.14) 
(2.8) and (2.9) imply that for A; = 0, 1, 

(1 - PNy'''idud.)''K^,o = (1 - e-^)dtKi-^j{u) + 4,ii^„,o (4.15) 
We first consider the case where A; = 0. Using (3.6), we get from (4.15) for all a; ^ / U J 

where A/uj is defined in (3.7). Setting 

Mr,,kiy^) ^ max (4.17) 

and, using that ^g'^^"(o"'' = ]E[T/|rf < rj], we obtain from (4.16) that for 1 — e~^(") < A/yj 
1 _ |1 _ e-"|MAr,o(^^)Mjv,i(0) < M^,o(^^) < 1 + |1 - e-"|MAr,o(u)Mjv,i(0) (4.18) 
But by Theorem 3.7 we have a uniform bound on Mjv,i(0), and this implies (4.13) for x ^ lUJ. 
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For k = I (4.15) gives 

(4.19) 

and the same arguments together with (4.13) for A; = show, for some c > and all 
< sR('u), |9(n)| < cc-^ryjj, that 

MnAu) < MN,m (1 + 0{\u\cnTjuj)) + |1 - e-"|M;v,i(«)Miv,i(0) (4.20) 
In particular, we conclude that on the same set, 

Mjv,i(n) = 0{MN,im = 0{cnTi^j) (4.21) 

Inserting this estimate into (4.19) (3.18) and (4.13) for A; = again gives for all < 
3?(tx),|9(u)| < ccjvTruj 

^^^ = (l + 0(HiVe^^--))^^^, x^/UJ (4.22) 

which yields (4.13) for A; = 1 and x ^ I U J. 

The remaining part, namely x € / U J, follows by first using (2.8), respectively (2.9), to 
express the quantities d'^Gj j in terms of d^G^j j with y / U J and then applying the result 
obtained before. <C> 

We now have all tools to establish a sharp relation between mean exit times and the 
principal eigenvalue A/ of P^. Set uj = —ln{l — Xj). We want to show that 

GZ,M) = 1 (4.23) 

Indeed, this follows from Lemma 4.1 with J = I U {m}, m G Mn, if we can show that 
\i < A/um- Now it is obvious by monotonicity that A/ < Xium- But if equality held, then 
by (2.6), lim^iuj G'm,/(^) — +°°! t)y continuity, it follows that there exists u < uj such that 
j(n) = 1, implying by Lemma 4.1 that 1 — e~" < A/ is an eigenvalue of P^, contradicting 
the fact that A/ is the smallest eigenvalue of P^f. We must conclude that Xj < A/um and 
that (4.23) holds. 

Theorem 4.4: Fix a proper nonempty subset IcMn- Let m G Mn\I be the unique local 
minimum satisfying Tj = Tm,i- Then 

Xi = il + OiTium/Ti)) E[Tr]-i (4.24) 
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In particular, 

Xi = RmTf^ (1 + 0(eiv|riv| + |riv|/(eivajv6jv))) (4.25) 



Proof: Using that for x > 0, > 1 x, for real and positive u, 

GZj{u) = E [e-^ ^^r-<r-] > P[C < + uE [CI.^<.;.] (4.26) 
Using this in (4.23), we immediately obtain the upper bound 

^ ^ P[rr<C] .427) 

Using now Lemma 3.8 to bound the right hand side, gives the upper bound of (4.24). The 
lower bound is of course already contained in Lemma 4.2. (} 

The a priori control of the Laplace transforms given in Lemma 4.3 can be used to control 
denominators in the renewal relation (2.10) which will be important for the construction of 
the solution of the equation appearing in (4.3). We are interested in the behavior of j 
near uj. 

Lemma. 4.5: Under the hypothesis of Theorem 4-4 there exists c > such that for all 

< K(U) < c/{cNTium) 

G^liu) - 1 = E [C]Ir„"J<rr] {U-UI + {U- Uj f 0{cNTium)) ^ ^ 

(4.28) 

= (1 + 0{eN))R-^ {U-Ui + {U- UifO{cNTium)) 



Proof: Performing a Taylor expansion at u = to second order of the Laplace transform 
on the left-hand side of (4.28) and recalling (4.23) we get 

G;^,,(n) - 1 = d^GZA^i) {{u - uj) -{u- ui)'ni{u)d^GZj{ujr') (4.29) 

where ^ 

ni{u) = [ sGZ /((I - s)ui + su)ds (4.30) 
Jo 

(4.29) then follows from Cauchy's inequality combined with (4.13) and (4.25) which shows, 
for c > small enough, C < cxd large enough, and all u considered in the Theorem, that 

(4.31) 

<C''cNTi^mduGZ,i{0) 
where we used Lemma 4.3. Using Lemma 3.8, the assertion of the lemma follows. <C> 
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5. Low lying eigenvalues 

In the present section we prove the main new result of this paper. Namely, we establish a 
precise relation between the low-lying part of the spectrum of the operator 1 — and the 
metastable exit times associated to the set Mn- Together with the results of Section 2, this 
allows us to give sharp estimates on the entire low-lying spectrum in terms of the transition 
probabilities between points in Mn and the invariant measure. 

As a matter of fact we will prove a somewhat more general result. Namely, instead of 
computing just the low-lying spectrum of 1 — P^, we will do so for any of the Dirichlet 
operators (1 — Pn)^ , with I C Mn (including the case / = 0). In the sequel we will fix 
IcMn with I ^ Mn- 

The strategy of our proof will be to show that to each of the points rrii e Mn\I corresponds 
exactly one eigenvalue Af of (1 — Pjv)^ and that this eigenvalue in turn is close to the principle 
eigenvalue of some Dirichlet operator (1 — Pjv)^S with / C Sj C A^jv- We will now show 
how to construct these sets Sj in such a way as to obtain an ordered sequence of eigenvalues. 

We set the first exclusion set Sq and the first effective depth Ti to be 

So = / and ri=rso (5.1) 

where Tk, KcMn, is defined in (1.9). If I / 0, let mi be the unique point in Mn\I such 
that 

Tm„i = Ti (5.2) 
If I = 0, let nil be the unique element oi Mn such that QAr(mi) = max^g_A4jv QNi'm)- 
For j = 2,... , jo, jo = \-M.n\I\, we define the corresponding quantities inductively by 

= T,j_2 U mj_i and Tj = Ty,-_^ (5.3) 
and rrij G MnX^j-i is determined by the equation 

Tiv(mj,Sj_i) = T,- (5.4) 

In order to avoid distinction as to whether or not j = jo, it will be convenient to set Tj^^i = 
b]^^. Note that this construction and hence all the sets Tij depend on N. An important fact is 
that the sequence Tj is decreasing. To see this, note that by construction and the assumption 
of genericity 

Ti = T^i,-Ei_-, > ejv^r„j+i,E;_i > ^nT^i+i,-Ei = ^nTi_^i (5-5) 
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The basic heuristic picture behind this construction can be summarized as follows. To each 
j = 1,... ,jo associate a rank one operator obtained by projecting the Dirichlet operator 
(1 — Pn)^^^^ onto the eigenspace corresponding to its principal eigenvalue As._^ ~ 
Note that our construction of Ej as an increasing sequence automatically guarantees that 
these eigenvalues will be in increasing order. The direct sum of these rank one operators 
acts approximately like (1 — Pn)^ on the eigenspace corresponding to the exponentially small 
part of its spectrum. Hence the difference between both operators can be treated as a small 
perturbation. 

Remark: We can now explain what the minimal non-degeneracy conditions are that are 
necessary for Theorem 1.3 to hold. Namely, what must be ensured is that the preceding 
construction of the sequence of sets is unique, and that the T^^ are by a diverging factor e^^ 
larger than all other T^^-Ej ■ 

We are now ready to formulate the main theorem of this section. Let Xj, j = 1, . . . , \Tn\I\, 
be the j-th eigenvalue of (1 — Pn)^ written in increasing order and counted with multiplic- 
ity and pick a corresponding eigenfunction such that {4>j)j is an orthonormal basis of 
I/c£2(r^^Q^). We then have 

Theorem 5.1: Set jo = \Mn\I\- There is c> such that the Dirichlet operator {1 — PnY 
has precisely jo simple eigenvalues in the interval [0, c6jv)|rjv|, i-e. 

a{il - PnY) n[0,cbN\T n\-^) = {Xi,. . . ,A,J (5.6) 

Define 7^ = oo and for j = 2, . . . ,jo 

= ^'^l^.Tmk,mj/Tj > ejf^ (5.7) 

Then 

A, = (1 + o{T-^ + r,+i/r,))) As,_, (5.8) 

where Xk, KgMn, is defined in (4-1) ■ 

Moreover, the eigenfunction satisfies for k = 1, . . . ,j — 1 

M^k) = ^jimj)0 {R^.T^^,^jTj) (5.9) 
Remark: Combining Theorem 5.1 with Theorem 4.4 and Theorem 3.5, we get immediately 
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Corollary 5.2: With the notation of Theorem 5.1, for j = 1, . . . ,jo that 



A, = (l + 0(^+T,+i/r,))E 
1 



= —Rmj (1 + C (|rjv|(ejv + l/(ajv6iveiv)))) 



(5.10) 



Note that Corollary 5.2 is a precise version of (ii) of Theorem 1.3. The estimate (5.9), 
together with the representation (4.4) and the estimates of the Laplace transforms in Lemma 
4.3, gives a precise control of the eigenfunctions and implies in particular (iv) of Theorem 
4.3. 

The strategy of the proof will be to seek, for each J = Sj , for a solution of the equation 
appearing in (4.3) with A near the principle eigenvalue of the associated Dirichlet operator 
(1 — P/v)^^-^. We then show that these eigenvalues are simple and that no other small 
eigenvalues occur. 

For the investigation of the structure of the equations written in (4.3) we have to take 
a closer look at the properties of the effective depths defined in (5.3). We introduce for all 
m e Mn\I the associated "metastable depth" with exclusion at / by 

rjv(m) = T^,MN{m), where Aljv(m) = I U {m' e Mn \ Qjv(m') > Qjv(m)} (5.11) 

Let us define for j = 2, . . . , _7o 

= min^Tmi,s,\m, (5.12) 

The following result relates our inductive definition to these geometrically more transparent 
objects and establishes some crucial properties: 

Lemma 5.3: Every effective depth is a metastable depth, more precisely for all j = 
1, . . . , jo ii follows 

Tj = rjv(m,)(l + 0(ejv|Miv|)) (5.13) 

For j = 2, . . . , jo we have 

> £j/Tj > e],\ (5.14) 
Moreover, for j,l = 1, . . . , jo, I < j, we have 

r^.EAm, =rsAm,(l + 0(eiv|Xiv|)) (5.15) 
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Proof: Fix I < j. It will be convenient to decompose E^- = U mi U S^, where = 
Ej\S;. Wc will use heavily the (almost) ultra-metric e(-, •) introduced in Section 2; for the 
purposes of the proof we can ignore the irrelevant errors in the ultra-metric inequalities (i.e. 
all equalities and inequalities relating the functions e in the course of the proof are understood 
up to error of at most ln3). Note that InT^j^j = e{x, J) — f{x), where f{x) = — InQAr(x). In 
particular, di = InT; = e(m/, — f{mi). As a first step we prove the following general 

fact that will be used several times: 

Lemma 5.4: Let m he such that e{m,mi) < e(m;,S/_i). Then f{m) > f{mi) + |lnejv|- 
Proof: Note that by ultra-metricity, 

e(m, = max (e(m, mj), e(m;, Ej_i)) = e(mj, S;_i) (5.16) 

But since for any m, 

e(m, Ei_i) - /(m) < di - |lne7v| = e(mi,Si_i) - f{mi) - \ \D.tN\ (5.17) 

which implies by (5.16) f{mi) < f{m) — |lneAr|."v' 

Let us now start by proving (5.14). The first inequality is trivial. We distinguish the cases 
where e(m;,E^) is larger or smaller than e(mi,Ej_i). 

(i) Let e(m/,E+) > e(m/,Ej_i). 

Since e{mi,'Sj\mi) = min (e(m/, E/_i), e(m;, E^)) , this implies that e{mi,Tij\mi) = 
e(m/,E/_i). 

Then, using (5.5) and genericity from Definition 1.2, 

e{mi,'Sj\mi) - f{mi) = e{mi, E/_i) - /{mi) = di > e(mj_i,E;_i) - f{nij-i) 

(5.18) 

> e(TOj_iEj_2) - f{mj-2) = dj-i > dj + \ Inejvl 
Obviously, this gives (5.14) in this case. 

(h) Let e{mi,T.^) < e(m/,E;_i). 

In this case there must exist nik G Ej" such that e{mi,T,j\mi) = e{mi,mk), and hence 
e{mk,mi) < e(mj,E/_i). Thus we can use Lemma 5.4 for m = m^. Together with the 
trivial inequality e{mk,mi) > e(mfc,Efc_i), it follows that 

e{mi,T,j\mi) - f{mi) = e{nik,rni) - f{mi) 

(5.19) 

> e(mfc,Efe_i) - /(mfe) -|- /(m;) - /(m^) > dk + |lnejv| > dj + \ \neN\ 
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This impUes (5.14) in that case and concludes the proof of this inequahty. 

We now turn to the proof of (5.15). We want to proof that the maximum over T^^s^Vrnj 
is reahzed for m = mi. Note first that it is clear that the maximum cannot be realized for 
m G Ttj\mi (since in that case T^^Y,j\mi ~ -'-)■ Thus fix m Sj. We distinguish the cases 
e{m,mi) less or larger than e(m, Sj\m;). 

(i) Assume e{m,mi) < e{m,Tij\mi). 

The ultra-metric property of e then implies that e{mi, T,j\mi) = e{m, Yij\mi), and hence, 
using the argument from above, /(m) > f{mi) + | Inejvl- Thus 

e{mi,'Ej\mi)-f{mi) = e{m,T,j\mi)-f{m)+f{m)-f{mi) > e{m,T,j\mi)- f{m) + \lneN\ 

(5.20) 

which excludes that in this case m may realize the maximum. We turn to the next case. 

(ii) Assume e{m,mi) > e{m,Tij\mi). 

We have to distinguish the two sub-cases like in the proof of (5.14). 

(11.1) e(m,,E+) >e(m,,Sj_i). 

Here we note simply that by (5.18) 

e{mi,T,j\mi) = e(mi,S;_i) - f{mi) = di > e(m, S^-i) - /(m) > e{m,T,j\mi) - f{m) 

(5.21) 

which implies that m cannot be the maximizer. 

(11.2) e(mi,S+) < e(mi,S;_i). 

This time we use (5.19) for some ruk £ S+ and so 

e{mi,T,j\mi) - f{mi) > dk > e(m,Sfc_i) - /(m) > e{m,T,j\mi) - f{m) (5.22) 

where in the last inequality we used that by assumption e(m, mi) > e{m, T,j\mi). Again 
(5.22) rules out m as maximizer, and since all cases are exhausted, we must conclude 
that (5.15) holds. 

It remains to show that (5.13) holds. Now the crucial observation is that by Lemma 5.4, 
MNirrij) n{m e Mn ■ e{mj,m) < e{mj,'Ej-i)} = (5.23) 
Thus, for all m G M.N{rnj), Tm^^m > rmj,Sj_i, which implies of course that 

Tmj,M{mj) > Tmj,Sj-i (5-24) 
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To show that the converse inequality also holds, it is obviously enough to show that the set 



(5.25) 



Assume the contrary, i.e. that for all m G M-{mj) Tmj,m > T^3,Sj_i- Now let m ^ / be 
such a point. Then also e{mj,m) > e(mj,Ej_i), and so by ultra-metricity e(m, Sj_i) = 
max (e(mj, m), e(mj, > e(mj,Sj_i). But, since /(m) < f{mj), it follows that 



(5.26) 



in contradiction with the defining property of rrij. Thus (5.25) must hold, and so TUj.,^jv(mj) ^ 
Tmj,i:j-i- This concludes the proof of the Lemma.<C> 

We now turn to the constructive part of the investigation of the low lying spectrum. Having 
in mind the heuristic picture described before Theorem 5.1 we are searching for solutions u 
of (4.3) for J = near us^.i = — log(l — As^._J. The procedure of finding n is as follows. 
The case i = 1 was studied in Theorem 3.5. For j = 2, . . . ,jo we consider the matrices 
Gj = Gi,^j defined in (4.2), i.e. 



ICj 

MY 



-9j 



l-G. 



and define 



\ 



■AG-^^i-^i, where = d\^g{l - GZ\^,:)i<i<j 



(5.27) 



(5.28) 



Equipped with the structure of the effective depths written in Lemma 5.3 and the control 
of Laplace transforms of transition times obtained in the previous chapter one simply can 
write a Neumann series for I — T>j{u)~^J\fj{u) for u near u-£._^ proving the invertibility of 
K,j{u). We then compute 



det = det ) = G,- det /C,- 



(5.29) 



where 



(5.30) 
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This follows by simply adding the column vector 

(-§)•) 

(which clearly is a linear combination of the first j — 1 columns of Qj) to the last column in Qj , 
and the fact that this operation leaves the determinant unchanged. Prom this representation 
we construct solutions Uj near itSj^i of (4.3). We begin with 

Lemma 5.5: For all j = 2, . . . , jo there are constants c > 0, C < oo such that for all 
C < CO and all 

CRm,£j^ <^{u)<cc-^^Tr^^, \'^{u)\<c/{cnTj+^) (5.31) 

the inverse of }Cj{u) exists. The l-th component of ICj{u)~^gj{u) restricted to the real axis is 
strictly monotone increasing and, uniformly in u, 

{ICj{u)-'gj{u))i = 0{mj\\u\-'Rm,T-lm, (1 = 1,... J -I) (5.32) 
Moreover, we obtain 

A = 1 - e"" G C7((l - Pjv)^) ^ Gj{u) = (5.33) 
where Gj{u) is defined in (5.30). 

Remark: Let us mention that the bound on SJ(u) in (5.31) is not optimal and chosen just 
for the sake of convenience. The optimal bounds with respect to our control can easily be 
derived but they are of no particular relevance for the following analysis. 

Proof: Fix j = 2, . . . , jo- Formally we obtain 

oo 

/C,(u)-i = (l-P(n)-X(u))"'P,(u)-i = Y,{V,{u)-^U,{u)YV,{u)-^ (5.34) 

s=0 

To use these formal calculations and to extract the decay estimate in (5.32) we must estimate 
the summands in (5.34). To do this we use a straightforward random walk representation for 
the matrix elements 

(5.35) 
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where to : mi ^ nik denotes a sequence to = (wq, . . . such that uoq = m,k, lo^^^ = ruk, 

ujt S U J) and uJt-i for all t = 1, . . . , |a;|. Assuming that the series in (5.34) 

converges, (5.35) gives the convenient representation 



aj:mi— >m, t=l 



where the sum is now over all walks of arbitrary length. We will now show that this sum 
over random walks does indeed converge under our hypothesis. 

By virtue of (5.15) we may apply (4.28) for m = mi and / = T,j\mi and conclude that 
there are c > and C < oo such that for all C < oo and all u G C satisfying (5.31) 

GZl^ji^) - 1 = (1 + 0{er,))R-] {u - us,\mj (l + - u^, vJOlcivTsJ) 

(5.37) 

= (1 + 0(ejv + 2c))uR-] 

where we used that u-£.\rni — cn^j- In addition, shrinking possibly c > in (5.31), (4.13) 
implies that for all A;, Z = 1, . . . ,j,k^l 

GZ,^,i^) = (1 + 0{\u\c^Tj+,)) GZ,^^{0) < 0(1)P[C^ < r^;] (5.38) 
Using these two bounds, (5.36) yields 

uj:mi — yrrij t=l 

To bound the product of probabilities, the following Lemma is useful: 

Lemma 5.6: Let a;o, wi, ^2,0;^ G such that coi 7^ coi+i, for all i and uq Uk- Then 



Proof: The proof is by induction over A;. For A; = 1 the claim is trivial. Assume that it for 
k = I. We will show that it holds for A; = Z + 1. Let s = max{0 <t<l\ujt = uo}. Note that 
by induction hypothesis and definition of s, 

i+i 

n ^[C-^ < ^S"] ^ < ^e;wa...vJ^[C.. < r^'M^^y—' (5-41) 

t=s+l 
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Now 



< T"s°\u;,+i\...vJ^['^<^!+i '^S^Aw,+i\...Vj+J 

■Wi+i ^ '(S3\wa+i\...\wi-i)Ua;j+iJ 



-^L'wj - 'SAu;«+i\...\wiJprT-t^i < r'^'l 



(5.42) 



Now the denominator on the right is, 
by (5.15). Thus, using the obvious bound 

s 

and once more that ujq G 'Ej\us+i\ ■ ■ ■ (5.42) inserted into (5.41) yields the claim for 

k = I + 1 which concludes the proof. 

Using Lemma 5.6 in (5.38) and the trivial bound R^^ < 1, we get 



(/c,(.x)-Vi(tx)).<P[C;<c;] ^ 

k=i ' ' ^ ' ' ^ 

< FW""' < r"^M — ^-^"^'1^1 — 

If C|Sj|f Itxl"-*- is say smaller than 1/2, the estimate (5.32) follows immediately. (5.33) then 
is a direct consequence of (4.3) and (5.29), since by (5.32) the determinant of ICj{u) cannot 
vanish in the domain of u-values considered. <0 

Remark: Defining 

Vj ^ diag(l-G-;_^^)i<K,„, Mi ^ Vi-gj,M^ and {f:f ^ (G-^,)i<fc<,o (5.46) 

where Gi,Mn defined in (4.2), a slight modification of the proof above shows that for c > 
small enough and all 3f?(u) < ch'^ such that 

^;K^JGZ,mM-A>{Vc)c-^' ^^T-]^^^^ (5.47) 
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one can write an absolutely convergent Neumann series for {l-VY^{u)Mi{u)) . Further- 
more, as a consequence of a random walk expansion similar to (5.45) we obtain the bound 

{Gi,MN{ur'fi{u))i = 0{aj'c],^T^^j) (5.48) 



This estimate is needed for the proof of Lemma 5.4. We are searching for solutions u near 
^^Sj_i of the equation appearing in (5.33). The case j = 1 is already treated in Theorem 3.5. 
Fix j = 2, . . . , Jo. We want to apply Lagrange's Theorem to this equation (see [WW]) which 
tells us the following: Fix a point a € C and an analytic function ^ defined on a domain 
containing the point a. Assume that there is a contour in the domain surrounding a such 
that on this contour the estimate |^'(C)| < \C ~ holds. Then the equation 

C = a + ^{0 (5.49) 

has a unique solution in the interior of the contour. Furthermore, the solution can be ex- 
panded in the form 

oo 

C = a + J^(n!)-^a^-^*(a)" (5.50) 

n=l 

We are in a position to prove 

Proposition 5.7: For j = 1, . . . ,jo there is a simple eigenvalue Xj = 1 — e~"^ < As^ such 
that (5.8), (5.10) hold if we replace Xj by Xj. Let (j)j he a corresponding eigenfunction. Then 
(5.9) holds if we replace (pj by (j)j. 

Proof: By means of Theorem 3.5 and (4.4) we may assume that j = 2,... ,jo. The equation 
in (5.33) can be written as 

G7„i,^{u)-l + ^j{0 = (5.51) 

where we have set ^ = '"^^^[''"s^^.J 

i-i 

^.(0 = J]G™;,^^.(n)(/C,(u)-V,(«))^ (5.52) 
1=1 

Fix constants c > 0, C < oo and let us denote by Uj the strip of all C € C such that 

Tj/£j < Sft(C) < cTj/Tj+,, |9(C)| < cTj/{Tj+,rNCN) (5.53) 

Putting Cs,_i = ^s,-iIE[rs/_J it follows Cs,_i = 1 + C(ejv) from (4.26) and (4.25) and we 
may apply (4.28) for c > small enough and all ( G Uj to obtain 

- 1 = E[T™^_J-^(l + 0{er,))R-] (C - Ci:,_, + (C - Ce,_ j'i?,(C)) (5-54) 
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where TZjiC) = E[r^.^_J-l7^s,_,(^x) is defined in (4.30). By (5.54) it follows that (5.51) is 
equivalent to 

C = Ce,_,+*,(C) (5.55) 

for some function '^j satisfying 

*i(C) = E[r^.^_J(l + OieN))R-]^j{0 + (C - Cs,_ J'7^,■(C) (5.56) 

Using (3.25) in combination with (5.4), it follows 

TZ.iO = OiT,+,/Tj) (5.57) 

Using (5.32) and the estimate (5.38), as well as (3.25), we see that for some c > 0, C < oo 
for all |C-Ce,_J < 1 

^jiO = E ^ (clTfT-l^T-l^,) < 0{c%Tr^) (5.58) 
1=1 

By means of (5.57) and (5.58) it follows for |C — Cej_i | < 1 

^j{0 = O{Tr^+Tj+,/Tj) (5.59) 

Since Tj > £j, by (5.14) and Definition 1.2, we may apply Lagrange's Theorem to (5.55) 
giving the existence of a solution Q = ■Uj]E[r^.^_J of (5.51) satisfying \Q — Cs^-J < 1- We 
rewrite (5.55) in the form 

= Ce,_, + 0{Tr^ + T,+i/T,) (5.60) 

By (5.33) Xj = 1 — e"^ defines an eigenvalue. Since from the invertibility of lCj{uj) it follows 
that the kernel of Qj{uj) is at most one-dimensional, (4.4) implies that Xj is simple. Using 
(4.24) and (4.25) for / = we derive from (5.60) that (5.10) and (5.8) hold, if we replace 

Xj by Xj. Moreover, using Uj < uy.^ from (4.4), we conclude that 

(^j(m/))i<Kj = 'Pj{mj)ICj{uj)-^gj{uj) (5.61) 

Hence from (5.32) and Uj = e'~'^^^uj]._-^ we obtain that (5.9) is satisfied if we replace (^j by 

4^j. 

Now it is very easy to finish the 



IE 
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Proof of 5.1: Proposition 5.7 tells iis that < Afc for A; = 1, . . . , jo- Assume now that 
there is /c = 2, . . . , jo such that Xk < Afc. Let /c = 2, . . . , j'o be minimal with this property. 
Since Afc_i = Afc_i is simple, we have Afc_i < Xk- Lemma 5.5 in combination with (5.30) now 
tells us that for j = 1, . . . ,jo some constants c > 0, C < oo and all Cc^Ej^ < u < cc'^^Tj~_^-^ 
the function Gj{u) is strictly monotone decreasing, i.e. has at most one zero. Hence from 
(5.33) for j = k — 1 and Gk-i{uk_i) = we deduce that uj, > cc~^^T^^. But since we already 
know that Uk < Cc^T^^ for some C, it then follows from (5.33) for j = k that Gk{uk) = 
implying the contradiction A^ = A^. 

Since Xj^ is simple, (5.33) for j = jo and Gjg{ujg) = implies Ajq+i > c6jv, where c 
denotes the constant appearing in (5.31). 

The remaining assertions of Theorem 5.1 then follow from Proposition 5.7. <^ 

6. The distribution function 

The objective of this chapter is to show how the structure of the low lying spectrum implies 
a precise control of the distribution function of the times r]", in cases where Theorem 3.5 
applies, i.e. IcMn, I,-Mn\I € MnXI, Tj = T^^,/. It is already known that 

the normalized distribution function converges weakly to the exponential distribution (see 
[BEGK] for the sharpest estimates beyond weak convergence in the most general case). 

The proof of these results proceeds by inverting the Laplace transforms G^iu), making 
use of the information about the analytic structure of these functions that is contained in 
the spectral decomposition of the low lying spectrum of (1 — Pn)^ obtained in the previous 
section. 

Let us denote by Cn the Laplace transform of the complementary distribution function, 

i.e. 

oc 

£Niu) = C^],{u) = J2 e'^'nrP > t] (5ft(«) < m), (6.1) 

where ui is defined in (4.26). The Perron-Frobcnius Theorem gives lim(l/t) logP[T™^ > 
t] = —uj. Hence the Laplace transform defined above is locally uniformly exponentially 
convergent. In order to obtain the continuation of £jv to the whole plane we perform a 
partial summation in the sum on the right-hand side of (6.1) and get 



GfUu)-l 
^ ^ e" - 1 



(6.2) 



38 



Section 6 



Invoking (2.8) a straightforward computation for A = 1 — e " shows that 

GIj{u) = {{1 - PnY - X)-\li^PNli)ix) (x^I), (6.3) 

Hence £jv is a meromorphic function with poles inu E {ui, ... , ii|rAr\7| }; where we recall the 
definition of the eigenvalues Xj = 1 — e""^ for j = 1, . . . , |riv\-?^| prior to Theorem 4.1. Since 
Cn is 27r-periodic in the imaginary direction, a short computation yields 

P[r7*i >t] = — : / e-^''CN{u)du. (6.4) 
27ri J_j7r 



Deforming the contour in (6.4) gives for Ujg < a < 'Ujo+i and C/q, = (0, a) x (— 7r,7r) 

e-*"£jv(u)dtx- 



-1 /"CK + iTT 

>t] = 7r^ e-'^CN{u)du - V e-*«^ res„, ^at, (6.5) 



where res„ Cn denotes the residue oi Cn at u. Here we have used that periodicity of jCn 
shows that the integrals over [a + Z7r,z7r] and [— ztt, a — iir] cancel and that the poles Uj, 
j = 1, . . . ,jo, are simple. 

Our main result can be formulated as follows: 

Theorem 6.1: Let jo = |A^jv\-f|. There is c> such that for some c> 0, 

F[tP^ >t] = -Y^ e"*"^ res„ . Cn + 6"*^=^^' {2Tri)-^ / e-*"£jv(^i)d^i, (6.6) 
j=i ' 

where uj = — ln(A7 + l) and Xj are the eigenvalues of {I — PnY that are estimated in Theorem 
5.1. Moreover, the residues satisfy 

res„, Cn = -1 + {Rm.CNT^/Ti) , res„, jCn = O {Rm.CNTj/T^) {j = 2,... , jo) 

(6.7) 

while the remainder integral on the right-hand side of (6.6) is bounded by 

{2TTi)-' [ e-'''/:N{u)du = 0{c]^^b]^^\rN\yn). (6.8) 

•J —in 



Remark: Recalling (3.25) and Theorem 4.1, one sees that Theorem 6.1 implies that the 
distribution of t^^ is to a remarkable precision a pure exponential. In particular, one has the 
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Corollary 6.2: Uniformly in t G E[r7*^]-iN 

F[Tp > tE[Tp]] = (1 + OiR^,CNT2/Ti)) e-<^+0{Rm,M/T,)) _ ^g ^^ 

We start with the computation of the residue of the Laplace transform at Ui . 
Lemma 6.3: 

res„, Cn = -1 + 0{R^,cnT2/T^). (6.10) 
Proof: Prom (4.23) for m = mi and the renewal relation (2.10) and (6.2) follows 

, _, gq.N 1 G7:^M^) 

res„, - Jim^ ^„ _ ^ GZ\M)-GZ\M~ e^^-ld^GZlM)' ^ ' 
Since ui = e'^'^^^ N-'^R^^T'^ , (4.13) for k = 0,l gives for some C < oo 

(u ) G™^ fO) 

- (1 + OiR^,CNT2/T2)) ^ '^Zl L ■ (6-12) 



Hence (6.10) follows from (6.11) in combination with (5.10) and (3.48). <0 

In general we cannot prove lower bounds for the higher residues for the reason described 
in the remark after Theorem 4.1. However, we can show that they are very small: 

Tiemma 6.4: 

res„. Cm = (Tj/T,)) (i = 2, . . . , jo). (6.13) 
Proof: For fixed j = 0, . . . ,jo we compute, using (6.2) and (6.3), 



u—^u 



1 u - Uj (ajcP/vij,(/)j)Q^ 

- e-«) 
<f>j{mi). 



AT = hm — — — (pjimi) 

6^^-1(1-6 "0-(i-e ") immN) 



(6.14) 

We can assume that (j)j{mj) = 1. We can express (pj(x), using the definition (4.4), Lemma 
4.3, and Theorem 5.1 in the form 

<l>j{x) =(1 + o(7))i^^,,E,(o) + j;o(r,/r^,^J(i + o(7))i^^„E,(o) ,^ ^ ^, 

1=1 (.o.lbj 
=(1 + 0(7))P[<. <rg^_J+0(7). 



40 



Section 6 



where 7 = max(T ^,Tj^i/Tj). Using Lemma 2.7, one sees easily that this imphes that 
for any e > 0, 

(I^Qn? > (1 + 0(e-^^))Qiv({x G r^v I |x - m,| < e/2}) > (1 - e)QN{A{mj)) (6.16) 
Prom (4.4) we conclude that, for J = Y!,j, 



k=l a^erjv 
j 



k=l »:erjv 
yei 



where we have used the symmetry of Pat. Applying (2.8) and (2.11) to the right-hand side 
of (6.17) we get 



k=l yel 



(6.18) 



, - V '7 > 



k=l 

Using that <pj{mj) = 1, we deduce from (5.9) and reversibility that 

Qiv(mfe)0,(mfe) = Qiv(m,-)0(i?-]r,/r^,,mj (6.19) 
Combining (6.19) with (5.38), (6.16), and, once more, (5.9) with k = 1, gives 



{m\Q.r^Mmi){T^ioPNT^i,cl>j)Q^ ^ 

k=l \ 



/ Tf 



(6.20) 



mi ,mj ruj ,1 



where we have used Lemma 5.6 for the sequences lo = {mj,mk, m) in the last equation. It is 
easy to verify that 

< rr. (6-21) 



Inserting (6.20) and (6.21) into (6.14), using uj = R^jTr'^ {1 + o{l)) and T^^jurm > Tj, we 
arrive at (6.13). 
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The last ingredient for the proof of Theorem 6.1 consists in estimating of the remainder 
integral in (6.6). This essentially boils down to 

Lemma 6.5: There is S > such that for all S'^R^^Tj^ < a < 6bN\^N\~^ o,nd all 
A = 1 — e~" on the circle |A — 1| = e~" we have 

Gl\{u) = 0{a-\-^'T^'). (6.22) 



Proof: Prom the strong Markov property (2.7) for J = 7 and L = M.n\I we obtain for 
3f?(u) < umn 

Kl,{u) = KIj^,{u) + Y,KT:MK'^umA^) e r^). (6.23) 

1=1 

Applying (1 — P^r — A)^ to both sides of the previous equation and evaluating the resulting 
equation at x = m/c, /e = 1, . . . , j'o, we conclude, as in (4.8), via (2.9) and (2.5) that 

= -Gl%,^{u) + Y^Gl\{u){5iu - GZ':,M»). (6.24) 
1=1 

Thus the vector 

^ (G-}(u))i<,<,, (6.25) 

solves the system of equations 

QiMn{u)^^ = fM. (6.26) 

where Qi^Mn^'^) fii"^^) defined in (4.2) and (5.46), respectively. In order to be able 
to apply (5.48) we claim that for some 5, c > 0, for all u = a + iv, v & [— vr, tt], and for all 
m G Mn\I 

We first observe that (2.2) shows that, for all SR(u') < umni 

<^N{m){GZ,MM - 1) = - J'iv)-^^^^ - A)Ki;^^,(n),<)^^(tx'))Q., (6.28) 

where we have extended the inner product to C'"'^ in the canonical way such that it is C-linear 
in the second argument. For |t; it 7r| < 7r/3 we simply get from (6.28), for u' = u and some 
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c> 0, using that a((l - Pjv)-^^\'")c(0, 1), 
|Q,v(m)Re(e-"(G-^^(n)-l))| 

= (((1 - P^)M.\m _ (1 + cos(r;)|)i^i,;^^(n),i^i,;^^(u) 

>(l + ce-«-l)(||<;^^(n)||Q,f 
>ce""Qjv(m). 

For |i; + 7r| > vr/S, — 7r| > 7r/3 and |f| > a, we derive from (6.28) for u' = u and some c > 



Qn 



(6.29) 



JAr(m)Im(e-"(G-^^(u)-l))|=|sin(x;)|e-'^(||Ki;)^^(«)||Q,) 

>Qjv(m)cae"". 



(6.30) 



In the remaining case, namely where \v\ < a, we use (6.28) for u' = u_x4N\m aiid obtain via 
(4.4), for I = MN\m, J = m, that 

\QN{m)e-^{GZ^j^Ju) - 1)1 = |A - A^,\„| Ki^^x. W'^Sa^.(^-M-w))q-I- (^-^l) 
Prom (4.13) it follows for some S > uniformly in a; G Fjv and \v\ < a 



(6.32) 



Since the minimum of the function |A — XmnXtu] is attained at A = 1 — e~", we conclude 
from (6.31) and (6.32) in combination with (4.4) for J = mi and (6.16) for some c > and 
all \v\ < a that 



lN{m)e-^GZ^^^{u) - 1)1 >c\X - A^„\™|(||k2>,, (^.M^v 

>c2Q;v(A(m))(l - e-"). 



(6.33) 



(6.33), (6.30) and (6.29) prove (6.27). Since by definition (5.3) and (5.14) it follows that 



(6.34) 



6jv is defined in Definition 1.1, combining (6.27) with (5.48) shows that the solution of (6.26) 
satisfies 

V'A(mi) = (V^a)i = O (a-^c-^^/Ti) . (6.35) 



Proof of Theorem 6.1: The proof of Theorem 6.1 now is reduced to the application 
of the Laplace inversion formula and estimation of the remainder integral. In view of (6.10) 
and (6.13) it remains to estimate the remainder integral on the right-hand side of (6.5). But 
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this is by means of (6.2) and (6.3) in combination with (6.22) for a = cb]\r\T ^ , < c < 5, 
fairly easy. (} 
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